Superconductor-Insulator Transition in a Disordered Electronic System 
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We study an electronic model of a 2D superconductor with onsito randomness using Quantum Monte Carlo 
simulations. The superfluid density is used to track the destruction of superconductivity in the ground state 
with increasing disorder. The non-superconducting state is identified as an insulator from the temperature 
dependence of its d.c. resistivity. The value of iTdc at the superconductor-insulator transition appears to be 
non-universal. 
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The problem of the effect of strong disorder on super- 
conductivity and of the resulting superconductor(SC)- 
insulator(I) transition in low dimensional systems has 
been studied experimentally for a number of years [f]. 
Theoretically, the problem is challenging because of the 
complicated interplay between interactions and disorder 
[2]. Within mean-field theory [3] superconductivity ap- 
pears to persist essentially all the way to the site localized 
limit due to an inadequate description of the disorder- 
induced fluctuations of the local order parameter. Much 
of the recent theoretical effort has thus focused on the 
dirty boson problem [4] which is expected to capture the 
essential physics of these fluctuations. The boson mod- 
els which describe the universal critical properties at the 
SC-I transition are also more amenable to analytical [4] 
and numerical [5.6] studies. However, if one is interested 
in characterizing the phases, and studying the question 
of a possible "universal" conductance at the SC-I transi- 
tion, one has to go back to a description in terms of the 
electronic degrees of freedom. 

As a first step in this direction we use Quantum Monte 
Carlo (QMC) simulations to study the simplest fermionic 
problem which can have superconducting, insulating, and 
(possibly) metallic phases. This is the attractive Hub- 
bard model with an onsite random potential 

-'^{l^ - '"i)ni<T - P\'^ni^ni^. (1) 
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We set t = 1 and measure all energies in units of t. Here 
Cio- is a fermion destruction operator at site i with spin 
(T, TiiCT = c|g.CiCT, and the chemical potential fixes the 
average density (n). The site energies are indepen- 
dent random variables with a uniform distribution over 
[— y, y]. The lattice sum {ij) is over near neighbor sites 
on a two dimensional square lattice. Note that this model 
focuses on the localization induced by the disorder; it 
does not, however, incorporate the disorder-dependence 
of the effective electron-electron interaction [2]. 



Our main results can be summarized as follows: 
f ) The transverse current-current correlation is used to 
compute the superfluid stiffness Dg. At low tempera- 
tures. _D.s decreases with increasing disorder T'. Beyond 
a critical the system becomes non-superconducting. 

2) For V > Vc the system is an insulator, as seen from 
the T-dependence of its dc resistivity p^^: dp,]r_/dT < 0. 
No evidence for a metallic phase is seen in 2D. A simple 
analytic continuation method, argued to be valid for dis- 
ordered systems, is used to estimate p^c for the first time 
using such QMC methods. 

3) The dc resistivity as a function of T for various disor- 
der strengths is also used to independently estimate the 
critical Vc for the SC-I transition. This estimate is in 
good agreement with the one obtained from D^- 

4) Pfir at V = Vr depends on interaction strength U and 
appears to be non-universal. 

Before turning to the results of the interacting, dis- 
ordered problem, we briefly summarize various limiting 
cases. In the noninteracting limit {\U\ = 0) one obtains 
the Anderson localization problem, which is known to be 
insulating for all V in 2D [2]. In the absence of disorder 
{V = 0) one obtains the attractive Hubbard model. Off 
half-filling, (n) 7^ f , the ground state has long range SC 
order and a finite temperature Kosterlitz-Thouless tran- 
sition [7]. QMC simulations have played an important 
role in the study of the non-random model, especially 
in its anomalous normal state behavior [8]. The same 
algorithms [9], which are still free of the fermion sign 
problem, are applied here to the model (1). 

We shall focus here mainly on various quantities that 
can be obtained from the current-current correlation 
function A [10]. The (paramagnetic piece of the) cur- 
rent operator is defined as 



J.(lr) = e^^ 



e-"\ (2) 



The Fourier transform of the impurity averaged A^x is 
then given by 
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A„(q;ic.„) = V / dT{j,{\,T)j^{G,G))e^^-'e-^--\ (3) 
J Jo 

where l^;,, = 2m: j (3, and {. . .) denotes a tliermal average 
at a temperature T = /3~^ for a given realization of disor- 
der and an average over an ensemble of such realizations. 
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FIG. 1(a). The kinetic energy, —Kx, and transverse and 
longitudinal current-current correlation functions, and 
A , are shown as a function of disorder V. tracks -A', 
as required by gauge invariance. The difference of A"'^ from 
A'a; signals the formation of a SC state at weak disorder (see 
Eq. (5)). 

Gauge invarianre recjuires that the longitudinal part of 
A satisfy the equality [10,11] 



lim, 
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A^^{q^,qy = 0;iuj„ = 0) = -K^, (4) 



where /i,, = (—ty^ (c!, - c, + cl c, , - )). the kinetic 
energy in the x direction, represents the diamagnetic part 
of the response. We have verified this ecpiality as a non- 
trivial check on our numerics; see Fig. 1(a). 

To study the destruction of SC with increasing disorder 
we look at the superfiluid stiffness obtained from the 
transverse current-current correlation function [10,11] 



A^ 
D, 



lim„ 



>o ^xx{qx = 0, Qy] iuj„ = 0) 



(5) 



7t[-K, - A^] 



Results [12] for A^, A^' and -A'^ are plotted in Fig. 1(a) 
for U = -4,T = 0.10, and (n) = 0.875 as a function of 
V. A"^ is estimated by using a linear extrapolation of the 
two smallest Qy values. In Fig. 1(1)) we plot U^, which 
decreases monotonically with disorder. There is a critical 
value Vc beyond which = and the system becomes 
non-superconducting. 

Another quantity of interest is the T = charge stiff- 
ness [10] D = 7r[— A'a; — lim^^o HeA^riq = O'.uJ + iO"'")], 
which is the strength of delta function d{Lo) in the optical 



conductivity. We find that D is indeed non-zero in the 
SC state and, within the accuracy of our numerics, we 
always find D ~ D^, as shown in Fig. 1(b). However, D 
is not useful in characterizing the non-superconducting 
state for V > Vc- In contrast to non-random systems, 
where D can be used to distinguish a metal (D ^ 0) 




FIG. 1(b). The superfluid density and charge stiffness 
D, as a function of disorder V . = D for all V. 

from an insulator (D = 0), for a disordered system D = 
for any non-superconducting state (since a dirty metal 
does not have a J-function in its T = optical conduc- 
tivity) . 

We must therefore turn to the dc conductivity to dis- 
tinguish a metal from an insulator. This is given by 
(Tdc = lim^_).o ImA(q = 0;uj)/uj, where A(q;w -|- 'iO+) = 
ReA(q;(j;) -|- ilmA{q;uj), and we drop the xx subscripts 
for simplicity. Using the fluctuation-dissipation theorem 
we obtain 



-t-cso 



duj exp(— wr) 
TT [1 — exp(—f3uj)] 



ImA^^(q;w), (6) 



valid for < t < /3. To obtain ImA from A(q; r), which is 
computed in the QMC, requires a numerical inversion of 
the Laplace transform. We will instead use a technique 
[8] which is valid for T <C fi, where fl is the scale on 
which ImA deviates from its low frequency asymptotic 
behavior (ImA ~ wudc)- Provided T <C Jl, Eq. (6) can 
be simplified to 



A,4q = 0;r = /3/2) = 7r(7dc//3' 



(7) 



which yields the dc conductivity. We note that this sim- 
plification may not be valid for non-random systems: e.g., 
for a Fermi liquid the scale ft ~ 1/Te_e ~ 7V(0)T^ so 
one can never satisfy T <C at low T. However, for the 
highly disordered state that we study, we expect the scale 
Jl to be set by the disorder V and to be T-independent, 
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so that Eq. (7) is valid. Wc will prosont below additional 
consistency checks of this approximation. 
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FIG. 2(a). pdc using Eq. (7) as function of temperature T 
for various disorder strengths V . At weak V, dp/dT > 0, but 
for large V dp/dT < (insulating). The point at T = 0.10 is 
obtained using Eq. (8) for V ~ V^. 



1 1 

U=-4.0 <n>=0.875 




o I 1 1 1 

2 4 6 

V 

FIG. 2(b). Pic as a function of disorder V for various tem- 
peratures T. Representative error bars are shown. 

In Fig. 2(a) we plot the dc resistivity p^c = l/fdc as a 
function of temperature for various degrees of disorder. 
We use units where = 7i = 1 so that the quantum 
of resistivity pq = h/(4e^) = it/2. For small disorder 
we see that p^c decreases with lowering T; that the un- 
derlying state is a SC can only really be seen from the 
superfluid density D^- With increasing disorder the T- 
dependence is altered qualitatively: for large V we see 
that dpdc/dT < 0, strongly suggestive of insulating be- 
havior (/5(jc = oo at T = 0). To see precisely where this 
transition takes place it is useful to replot the data of 



Fig. 2(a) as p^c versus the strength of the disorder V. 
with different curves corresponding to different tempera- 
tures. This is done in Fig. 2(b); from the crossing point 
of the various curves we estimate the critical disorder Vr 
separating the SC from the insulator. This crossing plot 
is not a consequence of any scaling ansatz, but simply 
follows from the monotonic behavior of p,]^ with T which 
we find in the QMC. These curves are remarkably simi- 
lar to those found in the experimental literature [1], and 
represent the first QMC calculations of the resistivity in 
a disordered, interacting fermi system. 

An independent estimate of Udc may be obtained as 
follows: We use the Matsubara conductivity 

a{iuj„) = [-K^ - Aj.j. (q, iw„)]/w„ (8) 

obtained from the simulation and fit it to a Drude form 
[15] ar)(/uj,-i) = (7dc/[l + l^')?!^]- which upon analytic con- 
tinuation yields cr£i{uj + iO'^) = CTdc/i^ — i^^T]. The simple 
Drude form is expected to be valid for a metallic sys- 
tem at low temperatures, where the uj and T-independent 
impurity scattering rate 1/r dominates any inelastic con- 
tribution due to the interactions. We thus use this pro- 
cedure to extract the low temperature = cr^^,^ for 
V Vc- The value obtained at T = 0.10 is shown 
in Fig. 2(a), and is in excellent agreement with the es- 
timates obtained from Eq. (7) at higher temperatures. 
cr(iujn) could be well fit by a Drude form only for V ~ Vc, 
suggesting the absence of an intermediate metallic phase 
between the SC and I in 2D. 

We estimate the critical Vc and the corresponding 
Pdc{Vc) from the raw data presented above, and sim- 
ilar data at other |L^| by 3 different methods. The 
first estimate of Vc is obtained from measurements of 
Ds at temperatures low enough that the pairing corre- 
lations are well formed across the entire lattice. Since 
~ (y - 14)^ with [4] ( = zv > 1, we expect to 
lie in the small tail in Fig. 1(b). The second estimate 
comes from the higher temperature crossing plots like 
Fig. 2(b), and the final estimate comes from the quality 
of the Drude fit, which works only for V ~ Vc, as de- 
scribed above. On the 8x8 systems studied, we find ex- 
cellent agreement between the three methods, as shown 
in Figs. 2(a) and 3(a). However, at this point we cannot 
rule out the effect of finite size corrections on our results 
[13]. 

We find evidence for only two phases in Fig. 3: a su- 
perfluid phase, with Ds > and p^i^ = at T = 0, 
and an insulating state with Dg = and p^c = oo at 
T = 0. We do not find any clear evidence in our nu- 
merics for a metallic phase (Dg = and = finite at 
T = 0). For \U\ = we know from the scaling theory 
of Anderson localization that there is no metallic phase 
in 2D [2] (although this wotdd be impossible to check for 
low disorder on the small systems under study here). In 
the opposite extreme |?7| — >• oo one obtains (short range) 
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ropnlsivo bosons in a random potential. There are gen- 
eral arguments [14] against the existence of a Bose metal 
with short range interactions at T = in any dimension. 
Thus, if there is a metallic phase in the 2D model we 
study, it can only exist in the intermediate \U\ regime, 
where the QMC should be most reliable. Note that the 
phase diagram for the same model (1) in 3D would cer- 
tainly be expected to have a metallic phase for small \U\. 
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FIG. 3. (a): Estimates for Vc at intermediate couplings 
|?7| = 3,4,6 from vanishing of (squares), and from p,-\c 
crossing (triangles), which have been offset for clarity. The 
full circle at the origin is the U = result = for 2D non- 
interacting electrons. We expect Vr (x t' /U at large couplings 
and ~ texp{-t/\U\) at small U. 

(b): Estimates for pdc(Vc) for \U\ = 3,4,6 obtained from the 
pdc crossing plots. 

While the results of Fig. 2(a,b) are in many ways rem- 
iniscent of experiments on the resistivity of thin films 
[1], we do not observe a "universal resistivity" [16] in 
the sense that pdc shows considerable variation with \U\. 
Recently, in a set of experiments on MoGe films, [l](d) 
a similar observation of sample dependent resistance was 
reported. In our model, \U\ directly controls the contri- 
bution of the unpaired electrons to the resistivity relative 
to the pairs. Our results for p^ir increasing with \U\ are 
consistent with a universal boson or pair resistivity and 
a fermionic contribution that increases with \U\. 
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